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Abstract 

We study the decay 7r+7r"e+e^, currently under analysis by the NA62 Col- 

laboration at CERN. In particular, we provide a detailed analysis of the Dalitz plot for 
the long-distance, 7*-niediated, contributions (Bremsstrahlung, direct emission and its 
interference). We also examine a set of asymmetries to isolate genuine short-distance 
effects. While we show that charge asymmetries are not required to test short distances, 
they provide the best environment for its detection. This constitutes by itself a strong 
motivation for NA62 to study K~ decays in the future. We therefore provide a detailed 
study of different charge asymmetries and the corresponding estimated signals. Whenever 
possible, we make contact with the related processes — tt^-k^^ and 
and discuss the advantages of A'+ — > Tr+Tr'^e+e^ over them. 



1 Introduction 



Kaon physics has played a crucial role in establishing the CKM flavor structure of the Stan- 
dard Model (see, for instance, [B [21 El llj ) • In the recent years, the experimental status of 
rare kaon decays has improved to the point that, at the time of writing, the NA62 physics 
program at CERN is aiming to detect eighty — Tr'^fiy events [5]. This constitues an 
extremely important test of the Standard Model (SM) and, as a result, a tool to constrain 
new physics [6]. The NA62 experiment aims to produce a total of 10^^ K~^^s, and hopefully 
maybe even 10^^ in the future. While many decays have been studied in the past, it is 
timely that studies be devoted to decays which have deserved less attention but that, 
given the experimental situation, can provide important information on CP violation, new 
physics and chiral tests. 

Traditionally, kaon decays have been divided in three categories: i) short-distance physics 
dominated decays, like the aforementioned — t- tt'^vu [3 [8]; ii) decays where short and 
long-distance contributions are comparable, like Kl — )■ TT^e~^e~ [9]; and iii) long-distance 
dominated decays, like radiative kaon decays [lOl EI]- However, even in this latter case, a 
particular charge, CP or angular asymmetry makes it possible to disentangle a small but still 
interesting short-distance contribution [61 |T2] . Some examples of this last category are the 
CP charge asymmetries studied by NA48/2 in K"^ — > 7r='=7r°7r° [T3j and K"^ — )■ 7r='=7r°7 [Hj. 

In particular, — )• 7r^7r*^7 requires a very accurate Dalitz plot analysis to disentangle 
the Bremsstrahlung contribution, accurately predicted from QED, from the direct (electric 
and magnetic) emission amplitudes [151 \TE\ [TTj . This Dalitz plot analysis is also useful to 
extract, directly from the NA62 measurement of the — t- tt^tt'^j charge asymmetry, the 
interference term proportional to the CP violating component of the electric emission piece, 
leading to a bound for the CP asymmetry, namely A^,p ^ < 1.5 • 10~^ |14j . 

In this paper we will study the decay — t- tt'^ ir^ e~ . This decay also belongs to cate- 
gory iii) above and is dominated by one-photon exchange, i.e., — )• 7r"'"7r''7* — )• 7r"'"7r*^e"'"e~. 
Its amplitude can be parametrized as 

Mld = J^^{k+,k^)^H^{pi,p2,q), (1.1) 

where pi are the pion momenta, ki the lepton momenta and is the invariant mass of the 
dilepton pair. is the hadronic vector, which, based on Lorentz and gauge invariance, can 
be written in terms of 3 form factors -Fi,2,3 as 

H^{Pl,P2, q) = F^pl^ + F2P^ + F3£'^'^"V.P2a'7/3, (1.2) 

where Fi^2 are dominated by the Bremsstrahlung (B) contribution, while -F3 contains the 
Magnetic (M) piece. These different contributions were studied in considerable detail in 
Ref. [18j. In particular, one finds the following rate for the branching ratios: 

5(K+ ^ 7r+7r°e+e") ~ B(E:+ ^ 7r+7r°e+e")ij + 5(i^+ ^ 7r+7r°e+e")M 

(330 ± 15) • 10"^ + (6.14 ± 1.3) • 10"^ (1.3) 

These numbers clarify how hard one must work to extract dynamical information from this 
decay: much like with — )• 7r^7r''7, one needs to fight the dominant Bremsstrahlung 
contribution. The numbers in Eq. ()1.3p were found by numerically integrating over the phase 
space using arbitrary kinematical variables |18] . One of the conclusions of the present work 
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is that an improvement can be obtained by studying the Dahtz plot with a specific set of 
kinematical variables, namely T* and E* (the kinetic energies of the charged pion and the 
photon, respectively, in the kaon rest frame) and the dilepton invariant mass q^. E* and T* 
were already introduced in Ref. [19] to study K'^ — )• ■k'^tt^'^ and thus adding comes as a 
natural extension. We will show that the set {E*,T*,q'^) allows to identify the most dense 
regions of the Dalitz plot for each dynamical contribution. 

Quite generically, as increases the relative weight of the Bremsstrahlung contribution 
gets reduced. Additionally, and most interestingly, for increasing q^ the different contributions 
tend to fill different regions of the {E*,T*) Dalitz plot. Thus, kinematical cuts in {E*,T*) 
can help discriminate the different dynamical contributions. Obviously, at high-g^ the phase 
space closes and thus experimental detection is afflicted with poor statistics. However, we 
will show that there exists a rather wide optimal window around ~ (50 MeV)^. 

The ability to gauge the values of is an evident advantage of — )• n'^'K^e'^e" over 

— )• 7r"''7r'^7, where one is restricted to q^ = 0, and thus exposed to a large Bremsstrahlung 
background. Even so, the Dalitz plot analysis of — t- tt^tt^'j has already been a useful 
tool to extract electric interference and direct emission amplitudes [14]. Besides improving 
on the accuracy of those determinations, — )• 7:^Tr^e~^e~ has a novel interference between 
Bremsstrahlung and direct magnetic emission. This provides a determination of the full F^, 
including the strong phases associated with final state interactions. 

In this article, we will also investigate in detail short-distance physics in this channel. 
Similarly to — ?• 7r"'"7r~e"'"e~[23l [Ml ESJ [26], short-distance effects can be obtained in 
— )• 7r~^ IT e~^e by looking into P-violation in the lepton pair. As a result, and opposed 
to — )• 7r^7r'^7, charge asymmetries are not needed for that purpose. However, in sharp 
contrast to Kl — )• 7r"'"7r~e"'"e~, in — )• vr^TT^e^e" there is no e— type contribution and, 
consequently, every CP-violating signal in the charged decay is a genuine short-distance 
effect. 

We will be interested in two different short-distance aspects, namely i) CP violation 
induced by — >• 7r"'"7r''7* and ii) short distances associated with effective dimension-6 
operators. For the latter, we will consider the contribution of SM physics as well as physics 
beyond the SM. In this last case, and in order to have a more predictive scenario, we will 
assume SM completions with Minimal Flavor Violation (MFV) [27]. We will show that, 
once MFV is assumed, rather precise statements can be made about short-distance signal 
detection. 

This paper is organized as follows: in Sect. [2] we will start with an analysis of the kine- 
matics. Sect. [3| is devoted to the dynamics of — )• 7r^7r''7*, which will then be used in 
Sect. [Hfor the study of the long-distance contributions to — )• vr+Tr^e+e". In particular, 
we will concentrate on the Dalitz plot of the photon-mediated P-conserving contribution in 
the {E*,T*) plane for different cuts in q'^. In Sect. [5] we will turn to the study of short 
distances with the evaluation of the SM contribution. Sect. [6] is devoted to the study of P 
and CP-violating signals in this decay within and beyond the SM. Conclusions are given in 
Sect, m 

2 Kinematical preliminaries 

In this Section we discuss the kinematics of K'^ — )■ TT^iT^e'^e~ in terms of the dynamical 
variables {E*,T*,q'^,9£,(j)), to be defined below. Since we are interested in a full Dalitz plot 
analysis, we will begin by providing the analytic expression for its phase space. 



2 



The invariant phase space for the four-body decay is defined as 

which can be triviahy rewritten as 

d^k+ Sk- 



(27r)32£;+ (27r)32£;_ 



5^^\k+ + k^-q), (2.2) 



where p-w = Pi + P2 and q = k+ + k- are the momenta of the dipion and dilepton pairs, 
respectively. Then one can obtain 

d<l> = ^(2vr)^ j ds^ j dseX^/\ml,pl, q^)<^^'^e, (2.3) 



where 



= / J2^,j2^/'^^P' ^P'-P-^ = (2^8^^'^'^^-"^'^'"^''^) /dcose., 

/ d'k+ d'k. + _ ) ^ ^1 /i _ M /d,/. /" dcos0,, (2.4) 

y (27r)32^+ (27r)32£;_ ^+ (27r)6 8 \/ q^ J J ^ ^ 

with A(a, 6, c) = a"^ + b'^ + — ab — ac — be. The angles above are defined as in Ref. [20j: if pi 
is the vr"*" momentum in the dipion CM system; k+ the momentum in the dilepton CM 
system; fi the direction of the dipion system as seen from the rest frame; and and k:}: 
the components of pi and k_|_ perpendicular to n, then 

" ■ Pi a ^ '^+ ^ Pi ■ ^+ 

cosC7^ = — — — ; cost'£ = — - — — ; cos0= — , . (z.o) 

IPiI |Pil|k+| 

Intuitively, 9i is the angle between the e"*" momentum and the dipion system as measured 
from the dilepton CM while 4> is the angle between the dipion and dilepton planes. 
The final result for the phase space therefore reads 



= \ 2 ^^^^^{'mj^,pl,q'^)^^^^{pl,ml+,mlo)dpldq^dcose^dcos 

(2.6) 

where the range of the kinematical variables is 

.2 ^ «2 ^ -.n^\'^ 



Am^ < q < {rriK — \Jp\) , 
{m„+ + m^of < pI < {niK - 2mif, 
< e^, de < TT, 

< (/> < 27r. (2.7) 

In order to disentangle most easily the different contributions, it is convenient to work 
with the set of variables {E*,T*,q'^,9i,(j)), where T* and E* are the kinetic energy of the 
charged pion and the photon, respectively, in the rest frame, and q^ is the dilepton 
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invariant mass. This set of variables is the natural extension of the variables introduced in 
Ref. [T9] to study — 71^11^^ decays. We can easily trade p^. and cos0,r for E* and T* 
using the relations: 



cos 9t^ 



where 



1 



Pit 



m\ - 2mKE* + g^, 
2iu-T*) 



rriK — E. 



2 V Pi 



dm 



(2i 



(2.9) 



and 5m = m^+— m^o is an isospin breaking parameter. In terms of the variables {E*,T*,q,9£, 
the phase space takes the simple form 



1 



2ii7r6 



/ 4m? ^ n ^ 
1 ^dEZdq^dT*d cos 



(2.10) 



and the limits of integration of Eq. (j2.7p . in terms of the set of variables {E*,T*,q'^,6£,(p), 



become 



u-^^Ef-q^ < T: < u+^^Ef-q^, 

mj^ - {m^+ + m^o)^ + q^ 



^ 2'mK 

2m£ < < rriK — ("i^+ + m^o), 

< 6*^ < vr, 

< (/> < 27r. 



(2.11) 



For later convenience we also list the following expressions: 



Q • P{1,2) = ^{mxE* - q^)C± =F rnK{u - T*), 
Pi-P2 = \{™^K - '^rnxE* + - {ml + m^,,)), 



mxC-l^e 



Q ■ P(i,2) = cos 0i/3l 



q^[mK{mK - 2EZ) + q'^][/3'^{Ef - q^) - 4(n - T*Y] sin sin 



/3[m;^(m;^ -2E*) + q^] ^ ^ [niKE* - q^){u - T*) 



-c±± 



AJEf - g2 



/3i sin 9i cos <J) q^[mK{mK - 2E*) + q^jE^^ - q^) - 4(n - r,*)^] 

2 Y _ ^2 



(2.12) 



which are needed to express the form factors and the kinematical weights (see next Section) 
in terms of {E*,T*,q^, 9^, (j)). 
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Figure 1: Photon-mediated — >■ 7r+7r°e+e decay with our kinematic conventions. The blob repre- 
sents the hadronic tensor H,, . 



3 Dynamics 

The amplitude for K+{P) 7r+(pi)7r°(p2)7*(9) ^ Tr+{pi)TT°{p2)e+{k+)e-{k-) can be parametrized 



as 



M 



LD 



■[u{k-)j''v{k+)]H^{pi,P2,q), 



(3.1) 



where if^ is the hadronic vector, which can be written in terms of 3 form factors -^1,2,3: 



(3.2) 



In terms of the momenta pi,2; q = k^ -\- k-. and Q = kj^ — k^, the squared amphtude can be 
cast as 



spins 



3 3 



(3.3) 



where 



Ti3 



q ■ Piq ■ Pj - Q ■ PiQ ■ Pj - q^Pi ■ Pj] = 1, 2), 

{Q-Pi)ef,uXpP'tp'2q^Q'; (^ = 1,2), 

Tss = 4mj [{ml - {pi ■ P2f)q^ - ml[{q ■ p2f + {q ■ pif) + 2pi ■ p2q ■ piq ■ P2] 
+ {Q • Pi? [{q • P2f - q^ml] + (Q • p2f [{q ■ pif - q^ml] 

+ 2Q ■ piQ ■ P2{q^Pi ■P2-q-Piq- P2)- (3.4) 

For phenomenofogical purposes, it is useful to split the form factors in terms of the (dominant) 
Bremsstrahlung and the direct emission contributions as Fi = F!> + F!> , where i = 1,2. 
According to Low's theorem |28] . the Bremsstrahlung piece can be written as 



A^(E:+ ^ vr+vrV)^ = 2e 



P-e 



+ 



Pi ■ e 



(P - g)2 - m\ {pi + qf - m\+ _ 



M{K+ ^+7r°) 



+ 



Pi 



{P-q)-\ {pi-q) + \ 



e^M{K+ ^ ^+^0), (3.5) 
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Figure 2: Diagrams corresponding to the Bremsstrahlung (upper row) and direct emission (lower 
row) contributions. Dotted and squared vertices represent 0{p^) and 0{p'^) weak operator insertions, 
respectively. Gauge invariance of the Bremsstrahlung contribution is ensured by Low 's theorem. 



from which one immediately concludes that 



iie{q-p2) 



{2q ■ pi + q^){2q ■ P - 



-2ie 



2q-P-q^ 



(3.6) 



The expressions for Ai(^K^ — )• vr+vr'^) and the direct emission piece (electric and magnetic) 
can be systematically computed using Chiral Perturbation Theory (ChPT) |18] . To leading 
order and neglecting isospin-violating effects they readj^ 



TT TT 



am 



(DE) 



Fl 



(DE) 



,{DE) 



;G27fn{mj^ - m^) - f.„6m { Gg + -G27 



2ieGse'^l 

In 

2ieG8e'^i 

In 

2eGHe'^i 



■P2MP+lq'M^^^+2q'L, 



(3.7) 



(3., 



^Knowing that BR{K+ -i> tt+tt") = 0.2066(8) and Ft = 5.32 • 10"^" MeV [8], one can easily work out that 
\M{K+ -> 7r+7r°)|2 = 3.28 • IQ-^^'GeV^ We will use this value in our analysis. 
^Regarding Gs and G27, in this work we adopt the conventions of Ref. [T]. 



6 



where 6q and 51 are strong phases associated with final state interactions (FSI) and we have 
defined the foUowing combinations of counterterms: 



<> 

(2) 



16 



Nu - iVi5 - N: 

- iVi5, 
iVi4 + 27Vi5 - 3(iVi6 

2 + 3(2a3 -aa; 



iV- 



17, 



N- 



17) 



1 



(3.9) 



Bremsstrahlung is a leading 0(jP) effect while direct emission appears as a subleading O(p^) 
effect. Notice that in the direct emission we have only included the counterterms: since 
the loop diagrams have no absorptive pieces their contribution can be absorbed in the Ni 
constants. The results above are analogous to the ones in K+ — 7r"'"7r'^7, except for the 



off-shell pieces in Fi^2- We note that above is also encountered in — )■ n^e'^e , while 

is a genuine contribution to — )■ vr+vr'^e+e^. These off-shell effects, as we will see 
later, turn out to be sizeable and make, even for rather small values of q'^, dramatic (and 
interesting) differences between — t- tt^tt^^ and K'^ — t- vr'^ 



- * 
vr 7 . 



4 Dalitz plot analysis 



where 



Here we will proceed in a similar way as it has already been done for — )■ tt ' 7r"7 
an accurate Dalitz plot analysis in the variables {E*,T*) was necessary to disentangle the 
Bremsstrahlung contribution, Tb, from the direct emission electric [Eqe) and magnetic 
(Mde) amplitudes. In that case one can write fW\ [T7] 



TT 



-vrS) 



dT*dE* 



dT*dE* 



1 + 2Re 





Ede 


2 


Mde 






cMk 


+ 







with 



and 



P 



e; 



Ede 



T* + m 



2 / 



1 



Mde 



E* 



rriK 



niK 



M ' 



rriK 



(4.1) 



(4.2) 



(4.3) 



Eq above is the vr'^-energy in the kaon rest frame and Aix is defined in Eq. (|3.7p . In compli- 
ance with Low's theorem, the Bremsstrahlung contribution can be factored out and one can 
extract the electric and magnetic parts from the p dependence. In this way NA48/2 obtained 
a determination of Ede from the interference term in Eq. ()4.ip together with jM/j^l |14j . 
Ede and Mde give, respectively, dynamical information on the non-anomalous and anoma- 
lous weak chiral lagrangian combinations Af^^ and which can be compared with 



existing theoretical predictions [291 [301 [3T] . Iii particular, it would be interesting to check 
independently the interference term, since NA48/2 finds opposite sign compared to the one 
expected on theoretical grounds. Moreover, — t- tt^tt^j has no reach over the sign of Mde, 
which has not been determined yet. 

Here we want to discuss — )• TT^ir^e^ e~ , complementing the kinematical variables used 
in Eq. (14. ip with the dilepton invariant mass, q^. As already emphasized in the Introduction, 
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(MeV) 


B [10-«1 

L J 


B/M 


B/E 


B/BE 


B/BM 




418.27 


71 


4405 


128 


208 


2 


307.96 


61 


3416 


111 


165 


4 


194.74 


48 


2320 


90 


129 


8 


109.60 


36 


1414 


71 


100 


15 


56.12 


26 


789 


56 


78 


35 


15.50 


16 


263 


41 


54 


55 


5.62 


12 


118 


38 


44 


85 


1.37 


9 


46 


49 


37 


100 


0.67 


8 


30 


71 


36 


120 


0.24 


8 


18 


458 


35 


140 


0.04 


9 


10 


-45 


37 


180 


0.003 


12 


5 


-19 


44 



Table 1: Branching ratios for the Bremsstrahlung and the relative weight of the rest of the contribu- 
tions for different cuts in q, starting at Qmin (fi-^st row) and ending at 180 MeV. In the last column 
we have also included the parity-odd magnetic- electric inteference term, to be discussed in Section,\d^ 



the — )• TT~^TT^e~^e~ decay rate is largely dominated by the Bremsstrahlung contribution 
{BR = 4.2 • 10~^), followed by the magnetic term (70 times smaller) and the electric inter- 
ference (120-130 times smaller). The pure electric term is much suppressed (40 to 100 times 
smaller than the magnetic). While sheer numbers show that the Bremsstrahlung dominates, 
this dominance is however not homogeneous in the kinematic variables. In Figure 3 we show 
the q dependence of the different contributions while Table 1 shows the integrated decay rates 
(the area under the curves) for different kinematic cuts in q. The general trend is that the 
Bremsstrahlung dominance gets reduced as q^ increases. However, even more information 
can be extracted if one looks at the distribution in the full {E*,T*, q'^,9e, (j)) hyperplane. The 
angular dependences are purely kinematic and can thus be integrated out. Eventually, the 
dynamically relevant objects we will focus on are the differential decay rates 

d^T d^Ts d^TE d^TM d^Tint . . 



dE*dT*dq^ dE*dT*dq^ dE*dT*dq^ dE*dT*dq^ dE*dT*dq^' 

where Tint collects the different interference contributions, BE, BM and EM. The decay rates 
above can be easily computed from 

d^Ti = ^ ^ IMldIU'^, (4.5) 



2mK 

spins 

where is the angular-integrated invariant phase space, whose explicit expression is given 
in Eq. ()2.10p . The matrix elements for the different contributions can be readily inferred 
from Eq. (j3.3p . We note that Tint above actually consists only of the BE term. This is so 
because the remaining (electric-magnetic) interferences are P-violating, i.e. odd in (f). Thus, 
both BM and EM terms cancel upon angular integration (c/. Section [6] for a full discussion 
of the electric-magnetic interferences.) 

For instance, with the definitions for the form factors and and the kine- 

matic weights of Eq. (13. 4p . and upon using Eq. (12.12p . one can easily show that the explicit 
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500 - 




0.00 0.05 0.10 0.15 0.20 

q{GeV) 



Figure 3: q dependence of the different contributions. The solid line represents the Bremsstrahlung. 
The dashed lines (from bigger dash to smaller dash) are 100 xM, lOOx BE and SOOxE, respectively. 

f 1 2) 

Error bars are omitted. The band around the electric contribution corresponds to varying Af^ ' in 
the range iA/'^"'' . The corresponding band on the electric interference term is completely negligible. 

expressions for the Bremsstrahlung and magnetic terms read (in the isospin limit): 



d^TB _ a'^\MK\'^{2mf + q'^) + \iq^ + Aag^ - 4m^c7 

dE*dT*dq^ ~ 487r3m3^g4 y ~ (^2 _ 2E*mKY{mK - 2{m^ + T* + E*))'^ ' 



dE*dT*dq^ 487r3/2g4 y q. 

where 



2 



Ai = 2mK{mK - ^E* - 2m^ - 2T*) - Ami, 

A2 = m^(16^f _ sE^rriK + {mx - 2{m^ + T*)){mK + 6(m^ + T*))) + IGmlniKE;, 
A3 = {mx - 2m-„){mK - 2E*) - 2E*mT,, 
A4 = iT*[mK -E;- 2m^ - T*], 
<T = A3 — mKimK — 2E*)X4, 

S^Sl-6l (4.7) 

Similarly, expressions for the remaining electric and electric interference terms can be found. 

In order to extract useful information from the previous expressions it is instructive to 
study the shape of each contribution in the (E*, T*) Dalitz plot while scanning over q'^. Thus, 
for any fixed q'^ we will have the corresponding Dalitz plot. This is a very convenient strategy 
to uncover how the different contributions populate the Dalitz plot and how this picture 
evolves dynamically. Throughout our analysis we will use the following set of parameter 
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values: 



mK = 493.677 MeV; 
a = (137.036)-^ 
Lg = 6.9 • 10"^ 

The values for A/"^"^ and have been taken from Ref. [14J. As mentioned in the Introduc- 
tion, the experimental value for turns out to be negative, in contrast to the theoretical 
estimates given in Refs. [29^ I30j . which yield a number comparable in magnitude but with 
positive sign. We want to point out that the holographic approach taken in Ref. [31] actually 
yields M^"^ = 0. This confirms that strong cancellations are a characteristic of M^'' and, 
therefore, that it is difficult to ascertain its value from hadronic models. The remaining weak 
counterterms Mj^ , in the absence of an experimental value, could in principle be estimated 
using hadronic models. However, a reliable estimate would require to evaluate the pion and 
kaon loops that we implicitly absorbed in the counterterms M^''^\ something that we will 

not do here. Instead, we will vary them in the range J^j^''^'^ G [~-^e^ t^e'^]- 

In Figure 4 we show snapshots of the different contributions at = (50 MeV)^. One can 
observe that the different contributions populate quite distinct regions of the Dalitz plot. At 
very low values of the BE piece is almost filling the same spot as the Bremsstrahlung, while 
the direct emission pieces stay in the center. This is precisely what one finds in — )• ir'^iT^'y 
and it is thus a cross-check that our results in the on-shell photon limit — )• are correct. 
However, even for modest values of q"^ the different contributions drift towards different 
edges of the Dalitz plot, whith the exception of the magnetic term, which remains centered. 
The picture shown in Figure 4 remains qualitatively valid up until values of close to 
the endpoint. As Table 1 illustrates, the rather distant but still looming Bremsstrahlung 
contribution diminishes its relative strength as increases. However, the overall signal to 
background ratio also drops sharply with g^, so for experimental detection large values are 
statistically disfavored. A fiducial regime can be found however at q^ = (50 MeV)^ or even 
smaller. This is possible because of the sharp evolution with q^ of the different contributions. 

We want to finish this Section by giving a quantitative estimate of the expected ex- 
perimental signals for all the contributions over the Bremsstrahlung, when one optimizes 
the information encoded in the {E* , T* ) Dalitz plot and concentrates on the regions where 
each contribution is most favored. Tables 2-4 show the expected ratios for BE, electric and 
magnetic contributions, respectively. The numbers quoted are the result of the integrated 
decay rates over a square centered at the peak of each contribution. Since the shape of the 
Dalitz plot changes with running q'^ (cf. the q^ dependence of E* and T* in Eq. (j2.1ip ). the 
squares are conventionally normalized such that their area is always 1 /8 times the full Dalitz 
plot areall The comparison is then done by integrating the same region for the remaining 
contributions. 

Table 2 shows the integrated Dalitz plot around the maximum of the electric interference 
term for different cuts in q'^. In order to see the improvement that can be achieved by 
placing cuts on E* and T* over a naive Dalitz plot integration, it is useful to compare with 
the results of Table 1. Notice that around the optimal cut = (50 MeV)^, the electric 

^We chose a square for simplicity. Any other shape (with the same area) would lead to similar results. 



= 137.293 MeV; 
U = 93 MeV; 

= -2.2164 • 10-3; 



5 = 5 



18' 



me = 0.511 MeV; 
Gs = 9.0 • 10-^ GeV" 
A/l? = 2.85 • 10-2; 



(4.8) 
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Qr (MeV) 


B/BE 




BE 




E/BE 


M/BE 


10 


701.90 


1.6 


• 10" 


-•n 


5.4-10-'^ 


0.07 


30 


165.34 


2.6 


• 10" 


-23 


2.3 • 10-3 


0.27 


50 


29.05 


1.0 


• 10- 


-23 


0.02 


1.23 


70 


40.75 


2.5 


• 10- 


-24 


0.01 


0.72 


100 


11.96 


6.9 


• 10" 


-25 


0.04 


0.60 


120 


7.07 


3.2 


• 10" 


-25 


0.07 


0.41 



Table 2: Branching ratios in a square centered at the maximum of the electric interference contribution 
(second column) for different cuts in q starting at qmin (fii'st row) and ending at 180 MeV. In the 
remaining columns we show the relative weight of the rest of the contributions. 



Qc (MeV) 


B/E 


BE/E 




E 


M/E 


10 


733.48 


49.26 


5.9 


10-24 


154.68 


30 


34.22 


4.32 


7.0 


10-25 


30.66 


50 


18.44 


0.71 


4.7 


10-25 


13.13 


70 


7.02 


-0.48 


3.1 


10-25 


5.27 


100 


3.57 


-1.15 


1.7 


10-25 


2.14 


120 


3.39 


-1.74 


1.1 


10-25 


1.37 



Table 3: Branching ratios in a square centered at the maximum of the Electric contribution (third 
column) for different cuts in q starting at qmin (fii'st row) and ending at 180 MeV. In the remaining 
columns we show the relative weight of the rest of the contributions. 



interference term is 30 times smaller than the Bremsstrahlung, which is comparable to what 
one finds in Table 1. However, electric and magnetic backgrounds have dropped significantly 
(from 25% to 2% for the electric and a factor 4 for the magnetic) . Comparison of Tables 1 and 
3 shows that the electric signal gets magnified by one order of magnitude with respect to the 
Bremsstrahlung. The electric intcrcfcrcnce background drops by roughly a factor 5, while, 
because of its proximity, no significant reduction of the magnetic background is achieved. 
Regarding the magnetic contribution, since the electric and interference backgrounds were 
already small, the important point is that it gets enhanced by a factor 2 with respect to the 
Bremsstrahlung. 

In summary, a much cleaner extraction of the interference and the magnetic terms is possi- 
ble: in both cases a significant enhancement of the signal is achieved (over the Bremsstrahlung) 
with small backgrounds from the remaining contributions. With this enhanced determination 
of the magnetic term one can even attempt an extraction of the electric piece: the signal is 
enhanced while the interference background gets reduced. 

As we mentioned above, we have neglected isospin breaking effects in our analysis. How- 
ever, notice that, since the Bremsstrahlung is the dominant contribution, isospin breaking 
effects there can constitute as sizeable an effect as the rest of the contributions. We have 
checked that isospin effects are sizeable only close to the Bremsstrahlung peak and therefore 
do not invalidate our previous conclusions. For completeness we list in the Appendix the 
explicit expression for the Bremsstrahlung contribution away from the isospin limit. 
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q,. (MeV) 


B/M 


BE/M 


E/M 




M 




10 


5.48 


0.34 


6.4 • 10-3 


9.1 


• 10" 


-22 


30 


5.21 


0.33 


9.6 • 10-3 


8.8 


• 10- 


-23 


50 


4.81 


0.30 


0.02 


2.5 


• 10- 


-23 


70 


4.46 


0.26 


0.03 


9.2 


• 10- 


-24 


100 


4.33 


0.10 


0.07 


2.2 


• 10- 


-24 


120 


5.30 


0.06 


0.10 


7.8 


• 10- 


-25 



Table 4: Branching ratios in a square centered at the maximum of the Magnetic contribution (fourth 
column) for different cuts in q starting at qmin (fii"st row) and ending at 180 MeV. In the remaining 
columns we show the relative weight of the rest of the contributions. 



5 SM short-distance contributions to K"^ ir^TT^e^e 

The effective Hamiltonian for AS* = 1 transitions reads [34^ I35j: 



|A5|=l^_G^ * 



y77(/^))Q77(/") + y7A{mw)Q7A{mw) + h.C, (5.1 



where Vij denote CKM matrix elements and 



Qjy = {msSRa^"'dL + mdSLaf"'dR)F^^, Qta = SLl'^di hi^lb^, (5-2) 

with L,R= 2^^^ ^^'^ ^i^^ ~ ~'''iltJ-^f ~ Qi^^fj.)- Below we will consider each operator in 
more detail. 



The first operator induces a CP- violating sd'j* vertex of the form |341 [35| [36] 

ie 

'W) 



nljj = V2GFj^V:,Vtd D'{xt)Q7y + h.C, (5.3) 



where 



It is a well-known fact [36j that the chiral realization of the hadronic matrix element in 
Eq. ()5.3p requires a dimensionally-suppressed operator at 0{p^). However, under renormal- 
ization this operator mixes with Q_ and its Wilson coefficient gets enhanced [37] This 
operator will be reconsidered later on in Section [6] when we discuss CP violation. 

Qta is induced by W-box and Z-penguin diagrams. In the presence of such contributions 
the matrix element for — )• tt'^ ir^ e~ takes the generic expression 

MsD = [u{k^h''-f5v{k+)]nf,{pi,p2,q), (5.5) 
where the hadronic tensor T-L^ is 

n>'{pi,P2,q) = + Fi^'pi^ + F|De^-"/^pi,p2a9/3. (5.6) 



'We note that even stronger enhancements can happen in BSM scenarios, e.g. supersymmetry 
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The interference between this axial leptonic current and the long-distance 7*-contribution 
discussed in previous Sections will give interesting P-violating effects that we will analyze in 
Sect. El 

In order to evaluate the contributions to the form factors F^^ above, one has to consider 
the bosonized version of the axial-vector effective operator 



^SD = C^^^=^S7/,(1 - 75)d 67^756 + h.c, 



(5.7) 



where si is a short-hand notation for the sine of the Cabibbo angle. The quantity ^ receives 
significant contributions from both top and charm quark loops and is given by [38^ [39] 



V*Vud 



where 



(5.8) 



(5.9) 



have been calculated including perturbative QCD corrections at the next-to-leading order, 
is of order 10~^ while is of order unity |40[ I41j. 

The chiral realization of the hadronic current in Eq. (15. 7p is, at the lowest order 



S7/.(1 - l^)d ^ -/^(L^)23 + ^E^.a.piL'L^Lf'U, 
where = iU^d^U . Thus, at very low energies one has 



T-LsD = iC- 



2^2 



(5.10) 



67/^756 + h.c. 

(5.11) 



Expanding out the U field above one easily obtains the following expression for the decay 
amplitude: 



M 



SD 



-a 



which immediately implies that 



P1-P2 



SD 



pSD _ 9 • 1-2 2 pSD 



Gpctsi 
2V2U 



(5.12) 

(5.13) 
(5.14) 



Notice that the three form factors are related by trivial numerical coefficients. This statement 
holds also when strong phases are taken into account: FSI in the previous form factors are 
described only by dl. Thus, in — t- TT^iT^e~^e~, 6q is exclusive of the Bremsstrahlung 
contribution while 6l is the only strong phase that affects short distances. 
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6 P- and CP-violation in — > Tx^Ti^e^e 

In the previous Sections we have concentrated on the (dominant) long-distance photon- 
mediated contribution to — )• 7r"*"7r'^e"^e~ as well as short-distance SM effects. In this 
Section we want to ascertain which asymmetries can be used to extract the different sources 
of short-distance physics. In order to carry out this short-distance analysis it will prove 
convenient to express the differential rate in the general form given for i^;4 decays [201 [21] : 



dE*dT*dq'^d cos 



= Ai+ A2 sin^ 6i + sin^ 6^ cos^ (/> -|- A4, sin 20i cos (/> 

+ ^5 sin 6i cos 4> + A% cos 0^ + A-j sin On sin 4> 

+ As sin 29i sin (j) + Aq s\v? 6e sin 2(f), (6.1) 



where di and (f) are two of the Cabibbo-Maksymowicz angular variables for K14 decays 
and Ai are dynamical functions that can be parameterized in terms of 3 dynamical variables, 
e.g., Ai{E*,T*,q). This factorized form makes it transparent which angular asymmetries are 
required to isolate the different Ai. 

One can easily show that the Bremsstrahlung, direct emission and electric interference 
terms contribute to A1-4. In contrast, ^8,9 receive contributions from the electric-magnetic 
interference terms (BM and EM) and therefore capture long-distance induced P-violating 
terms. -4.5^6,7 also P-violating terms but generated through the interference of Qja with 
long distances. 

In this Section we will extract information about ^5-9 through a set of angular asym- 
metries. We will thus discuss P violation both at long and short distances, but also the 
extraction of CP- violating signals from — )• e"^ e~ . 

6.1 P violation at long distances 

While doing the Dalitz plot analysis in Section [3] we mentioned that the magnetic interference 
term cancels after a full angular integration. This is due to the intrinsic parity-odd nature of 
the magnetic term. Actually, integrating Eq. (16. ip with respect to Oi, the differential cross 
section looks like: 

— = Ii cos^ 0+2^2 sin^ 0+2^3 sin cf) cos <f> +I4 sin <f> +I5 cos (f), (6.2) 
dcp 

where Xi 2 are P-conserving and include the B, E, M and BE contributions while X3 4 5 are 
P-violating. In this Section we will be interested in X3, which consists entirely of the BM 
contribution and thus probes P violation at long distance^. In contrast, ^4^5 are related to 
^5,6,7 and therefore test P violation at short distances (c/. the next Section). 

One way to isolate X3 is to define an appropriate asymmetry in (p. Since this term is 
asymmetric with respect to the shift (p ^ ir — (p, one defines the following piece-wise 
angular integration 

r2TV r /-tt/S /-tt /■37r/2 p'^n ' 

(6.3) 



/>27r 


r7v/2 pn p3n/2 /'27r 


/ = 




Jo 


Jo Jit/2 Jit Jztt/2 



^The EM contribution, also inside Z3, is negligible. 
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(MeV) 


B/BM 


BE/BM 


E/BM 


M/BM 


BM 


10 


593.03 


1.57 


2.0 • 10-3 


0.31 


2.6 • 10-2^ 


30 


150.75 


1.56 


6.6 • 10-3 


0.88 


2.9 • 10-23 


50 


85.01 


1.45 


1.5 • 10-2 


1.46 


7.6 • 10-24 


70 


58.77 


1.20 


3.4 • 10-2 


2.02 


2.5 • 10-24 


100 


40.29 


0.55 


0.12 


2.78 


5.4 • 10-25 


120 


33.90 


-0.13 


0.26 


3.17 


1.9 • 10-25 



Table 5: Branching ratios in a square centered at the maximum of the magnetic interference contri- 
bution (fifth column) for different cuts in q starting at qmin (fi-Tst row) and ending at 180 MeV. In the 
remaining columns we show the relative weight of the rest of the contributions. 



then the P- violating observable 



27r 



dV 



4'^ = 'WJr^ ^ (6-4) 



f^'^ dV 
Jo d(f) 



naturally selects the magnetic interference piece. Indeed, the contributions for Xi, X2, X4 and 
X5 identically cancel above and we are left with 

d^TBM ^ _^GsMkM^ L _ 4m[\ ^'^ g4 ^ (2A3 - \A)q^ + a ^.^^ 

dE;dT*dq^ 247r4/^mK(?2 V q^J ^ E*^^ - q^{mK - 2{E* + + T*)) ' 

(6.5) 

where Mr, J^e^ and J\f^^ are defined in Eqs. (j3.7p and (j3.9p . while A3, A4, a and 5 are 
given in Eqs. (14. 7p . Incidentally, we want to note that the structure entering the numerator 
coincides with the one found in the magnetic term (see Eq. (j4.6|) ). 

The observable defined in Eq. (j6.4p was already used to extract the magnetic interference 
term in Kl — >■ T:^'K~e^e~ [351 [33]. However, in that case the 0{p'^) contribution is strongly 
suppressed and 0{p^) terms dominate. In — )• -K'^ir^e^e" , in contrast, the 0{p^) is 
already sizeable and counterterms in are expected to be small. A measurement of 

A^p^ can confirm the sign of the magnetic interference and thus the convergence of the chiral 
expansion. Notice that the previous contribution can also be used to extract the strong phase 
5. 

In Figure 5 we show the Dalitz plot at the representative g = 50 MeV cut. We note 
that the shape of the Dalitz plot for this contribution is largely independent of q^ (this 
is to be expected because neither the magnetic nor the Bremsstrahlung pieces vary much). 
Table 5 shows the comparison of the expected signal against the (Bremsstrahlung-dominated) 
background. As discussed in previous sections (c/. also Table 1), g = 50 MeV is probably the 
optimal cut: lower values of (f' are largely dominated by the Bremsstrahlung while higher 
values are afflicted with poor statistics. Therefore, naively (without better cuts), the precision 
required to extract the BM interference is of the same order needed to isolate the magnetic 
term in — >• ■k'^ ir^ e'^ e~ . 
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6.2 P violation at short distances 



In this section we would like to isolate P-violating effects that come entirely from short 
distances. This can be done by considering P violation not in the hadronic tensor but in the 
dilepton pair. This way one ensures that there is no contamination from the long-distance 
QED background [23]. We have already noted in Section [5] that axial leptonic vertices, 
which are generated by VF-box and Z-penguin diagrams, induce P-violating effects through 
interference with vector leptonic currents. In this Section we will use this fact to study P 
violation both within and beyond the SM. 

Similar to Section [5] we will consider the generic expression: 



where T-L^ takes the form. 



MsD = [u{k.)l^-i^v{k+)\H^{pi,p2,q), (6.6) 



= Tipt + J-2P^ + Tze^'''''^PiuP2aqp. (6.7) 



The form factors J-i contain short distance information, namely the W and Z one-loop con- 
tributions evaluated in Section [5] as well as potential contributions from physics beyond the 
SM. With the normalization we have chosen, = Ff^ + Tf^^'^. 

The dominant P-violating piece will arise from the interference between the long-distance 
(photon-mediated) matrix element Mld and the short-distance AisD with an axial dilepton 
current. The resulting amplitude reads 

MldM*sd + MsdMId = ^Im [H^^n"*] e^.xpklk'L = %Y. " ^3^11 (6-8) 

where Fi are the form factors defined in Section [3] and Gij are given by 
G\2 = f^f,u\pPiP2Q^q'', 

Gi3 = -[pi- Q{q^Pi ■P2-P1- qP2 ■q)+P2- Q{{pi ■ qf - ml+q^)] , 

G23 = [P2 ■ Q{q^P2 ■P1-P2- qpi ■q)+Pi- Q{{P2 ■ qf - "i^og^)] . (6.9) 

As we remarked in the Introduction, the form factors are functions only of the variables E*, 
T* and q"^, while the angular dependence is entirely contained in the kinematical weights Gij. 
Therefore, even without knowing J^i explicitly, we can readily isolate the angular dependence 
of the vector-axial interference term as 



(fT 

= A5 sin 0i cos (j) + Aq cos Oi + Ay sin 9i sin cj), (6.10) 



dcos 

where A^^q receive contributions from the G13 and G23 terms in Eq. (|6.8p . while comes 
entirely from the term proportional to Gi2- Notice from the Gij above that the vector-axial 
interference is linear in Q. An appropriate asymmetry to be sensitive to such contributions is 
therefore the dilepton momenta asymmetry k± ^ —kzp. In terms of angles this corresponds 
to {6e, (/))—)• (vr — Oi, n + (p), and thus a possible observable is 

1 r^/^ (p-p rO |.37r/2 ^ ^2p 

d COS Of / dS—— / dcos6f 



AS) ^ Jo Jo dcpd cos 6 e J-i Jj, dc/xicos 6e , . 

P ~ rl /■1-/2 £y /"O /■3'r/2 ^2p ■ l"-J^J-J 

d COS 6 1 I d(f)-— —+ I dcosOg, 



Jo 



d(j)d cos 9 1 J_i d(j)d cos 01 
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Ap selects, by construction, short-distance information: BSM or weak interaction couplings 
of the leptons. For definiteness, and in order to be conservative, we will assume Minimal 
Flavor Violation (MFV) \27\ for BSM interactions. This means that no new operators other 
than the SM ones will appear. Consequently, the form factors F are constant and satisfy the 
relations of Eq. ()5.13p . Notice also that MFV implies the existence of only one weak phase 
and one strong phase. Taking all this into account, the asymmetry of Eq. (|6.11|) takes the 
form 



where F(^) = ' + F^ ' , Mde is given in Eq. (fOj) . 5 is the FSI and the brackets stand 
for phase space integration. Use has been made of Eqs. (I5.14p to express the short distance 
form factors in terms of only Fi and J-3. This way the magnetic or electric character of the 
form factors is transparent. 

Some comments are in order: 

• Since is not a charge eigenstate, Ap contains both CP conserving and CP violat- 
ing contributions. The latter (proportional to G12) are suppressed and therefore not 
included above (see, however, next Section). 

• The two terms proportional to Fi are competitive: Fi and Mp)E have a strong phase 
difference of 90 degrees and therefore there is no strong phase suppression in the second 

term. This compensates for the suppression of |-/Vf£)£;| over IF^^^^I. However, it is the 

is) 

term proportional to J-3 which eventually dominates Ap . 

In Figure 6 we show the Dalitz plot for each contribution, always assuming MFV, i.e., constant 
Fi- Only the shape should be considered meaningful: without specific model predictions for 
the form factors Fi, the plots have been normalized to arbitrary units. However, they point 
out the most favored regions to search for each source of P violation. As a matter of fact, 
the second term in Eq. ()6.12p is positive definite, while the first and third terms are not. The 
plots can be qualitatively understood by noticing that the first and third terms are akin to 
Bremsstrahlung-electric interference (modulated by the kinematical factors Gjj), while the 
second one is a magnetic-electric interference (again modulated by Gij). 

6.3 CP violation 

In the previous Section we studied direct contributions to short-distance physics through 
genuine short-distance operators. One of the conclusions we reached is that new physics can 
be detected without resorting to charge asymmetries. This is particularly welcomed since 
no K~ detection is needed to explore short distances. However, we will show that when the 
charge asymmetry can be built much more short-distance information can be extracted. 

In order to have a CP-violating charge asymmetry in — )• vr^vr^e^e" one has to choose 
a P-conserving observable. In contrast to the neutral mode Kl — )• ■n'^ ■n~ e'^ e~ , here no long 
distance (e-like) contributions are present and thus charge asymmetry is directly sensitive 
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to short distance effects (electroweak and/or BSM phases). CP-violating phases might show 
up either associated to new physics operators (c/. last Section) or induced through — )■ 
TT'^TT^'j*. Therefore, there are two possibilites: i) terms proportional to (P-conservation in 
both the leptonic and hadronic pieces) or ii) terms proportional to ec^ (P violation in both 
the leptonic and hadronic pieces). In this Section we will consider both of them. 

The simplest observable one can build is the fully integrated charge asymmetry, which 
yields the result 

r(E:+ ^ 7r+7r°e+e-) - r(i^- ^ 7r-7r°e+e-) eV,^,,^ A . ^ 
^CP = ZT7^ + + , _ + _s ~ Eb\\Ede\ )smdsm<^E, 

(6.13) 

where we have defined Eb = A4{K~^ — 7r~^7r^'j*)B {cf. Eq. p.Sp ). This interference already 
appeared in — )• tt^tt^j, with the lower bound ^ < 1.5-10"^ [H]. We want to remark 

that, exactly as in — )• 7r^7r'^7, the appearance of the SM weak phases $e,m occurs at 
0{p^) in the chiral expansion of the direct emission form factors |36l I37j . Therefore, in 
MFV, E£)E is not the one given in Eq. (j4.2p but is related to the chromoelectric operator 
Q^-^ discussed in Section [5l However, we note that 0{p'^) contributions can be generated in 
non-MFV scenarios [43 1 i44j . 

The main advantage of — )• T:^T:^e^e~ is that one can also build a CP-violating 
observable not afflicted by strong phase suppression. This can be achieved through a charge 
asymmetry with the angular asymmetry discussed in Section [6.11 

^CP = '^'^ ^(\Eb\\Mde\) cos5sin$M, (6.14) 



(K++K-) , , B \ / 



-asym 



where the magnetic amplitude Mde is, again, a chirally-suppressed 0{p ) contribution [36| 
I37j . Notice that this is a genuine CP-violating contribution of — )• TT^TT^e~^e~ , which 
comes from the fact that the virtual photon can generate different helicity amplitudes. We 
want to remark that this term might give a CP-violating signal 2 orders of magnitude larger 
than the one in K+ —7- tt'^tt^'j, i.e. Aqp (i) it avoids strong phase suppression and (ii) 
a better detection can be achieved focussing on the central part of the Dalitz ploti where, 
for moderate values of q^, there are strong enhancements of the signal over background ratio 
(see Table 4). 

The observables Acp and A^p considered before select CP violation in 7-mediated — t- 
TT^T^^e^ e~ . One can also access CP violation in direct short distance operators to — t- 
TT^TT^e~^e~ by defining an alternative ^-asymmetry, namely 



2lT 







tt/2 f-K />37r/2 f2TT 

+ 



Jtt/2 Jtt J3tt/2 



(6.15) 



We are here assuming that the magnetic contribution induced by Qt^. does not differ kinematically from 
the 0{p'^) one considered in Section [S] 
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which naturahy leads to the observabl^ 



^CP = %.,r t^{Gi2\T,\\Fb\ ) cos<5sin$. (6.16) 



Notice that A'^p actually selects the CP-violating piece of Ap^ defined in the previous Section. 
There it was competing with CP-conserving terms and, since it was a subdominant effect, it 
was neglected in Eq. (j6.8p . 

Due to the structure of the NA62 detector, only the forward direction in cp can be accessed. 
This suggests to define the forward asymmetry: 



+ 



37r/2 



Notice that, in view of the observables defined in Eqs. (j6.14p and (j6.16p . they both contribute. 
Therefore, this forward asymmetry is a way to optimize the detection of direct CP-violating 
phases. 

It is instructive to compare the situation with Kl — )• vr^vr^e^e". There CP observables 
are dominated by indirect CP violation, which overshadows any direct CP contribution. One 
way to isolate direct CP violation was put forward in Ref . [25] , where the analogue of Ap 
was considered. The advantage of — t- TT~^TT^e~^e~ is that direct CP violation is free from 
indirect CP violation backgrounds. This, added to the fact that NA62 will collect decays 
with large statistics, makes the study of this channel appealing. 



7 Conclusions 

In this work, we have studied different aspects of the decays — )• 7r^7r*^e"^e^. At long 
distances, we have performed a thorough Dalitz plot study. In the past, cuts in the leptonic 
invariant mass were shown to be useful to extract information about the different contri- 
butions of this decay. Here, we have shown that additional cuts on the variables {E*,T*) can 
provide even cleaner signals for a precise extraction of all contributions by concentrating on 
different corners of the phase space. 

In particular, interesting information about the direct emission (electric and magnetic) 
contributions can be extracted through their inteference with the dominant Bremsstrahlung 
contribution. The electric interference was already determined from — )■ tt~^tt^^ by the 
NA48/2 collaboration, leading to a discrepancy in sign with the theoretical estimates. — )• 
7r"'"7r''e~''e~ offers an alternative measurement with the potential to sort out the discrepancy. 
The magnetic interference term is however a genuine feature of — )• TT~^TT^e~^e~ that will 
allow for a measurement of the sign of the magnetic term, something relevant for chiral tests 
of the anomalous weak chiral Lagrangian. Additionally, it provides a clean determination 
of the FSI phases in — )• vr^vr^. In that sense, K+ — )• TT~^7r^e~^e can be considered a 
laboratory for chiral tests. 

^The terms proportional to G13 and G23 can in principle be accessed with a different angular asymmetry, 
but they are suppressed by strong phases. 
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Short distances within and beyond the Standard Model can be probed in this decay 
channel without resorting to charge asymmetries: observables with P violation in the dilepton 
pair coupling are already a signature of short-distance physics. Thus, even at this stage of the 
NA62 experiment (with only K~^), short-distance analyses can be performed. Regarding the 
charge asymmetry, an interesting feature of — )• TT^Tr^e~^e~ over its neutral counterpart 
Kl — 7" TT~^TT~e~^e~ is that short-distance observables are not contaminated by indirect CP 
violation. We have also explored the potential of charge asymmetries to uncover new physics, 
should NA62 eventually start data recollection on K~ . With the help of the Dalitz plots and 
assuming Minimal Flavor Violation, a very conservative BSM scenario, one can show which 
are the optimal regions to look for signals of CP violation. 
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A Appendix 

For completeness we here list the differential decay rate for the Bremsstrahlung contribution 
of Eq. ()4.6p with isospin breaking corrections taken into account: 

a^\MK\'^{2mf + q^) I 4mf + X^^q^ + A|g^ + 

ASTT^mj^q^iq'^ -2E*mK)^\l q^ [(m^ - 2(m^+ + T* +£;))- x(m^ - 2m^)]2 ' 

(A.l) 

where we have defined 

= Ai - 2mj^x - mj^X^, 
A^ = A2 + 2mlcX {rriK - 2mK{2m^ + 3T*) + 4m^(m^ + T*)) , 

+ it^kX^ {AE*mK - 3m\ + 4mKm.T, + 4m^) , 
A^ = -4mla - Sm]^\^x{E; - rriK + 2(m, + T*)) - Am%x^{E; - mK + 2m^f , 

and the isospin-breaking parameter 

— m^o 

X = • 

rriK 
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Figure 4: On the left panels: Dalitz plot in the {E*,T*) plane at = (50 MeV)'^ for the different con- 
tributions: B, BE, E and M (from top to bottom). On the right panels we show their two-dimensional 
density projections. 
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Figure 5: On the left panel: Dalitz plot in the {E*,T*) plane at = (50 MeV}^ for the P-violating 
BM contribution, obtained after a piece-wise angular integration. On the right panel we show the 
corresponding two-dimensional density projection. 
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Figure 6: Dalitz plot for the contributions to Ap (in arbitrary units) at = (50 MeV)^. The 
different terms are geometrically placed as they appear in Eq. i6.12\) . 
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